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ABSTRACT 

Explosions of sub-Chandrasekhar-mass white dwarfs are one alternative to the standard Chandrasekhar-mass 
model of Type la supernovae. They are interesting since binary systems with sub-Chandrasekhar-mass primary 
white dwarfs should be common and this scenario would suggest a simple physical parameter which determines 
the explosion brightness, namely the mass of the exploding white dwarf. Here we perform one-dimensional 
hydrodynamical simulations, associated post-processing nucleosynthesis and multi-wavelength radiation trans- 
port calculations for pure detonations of carbon-oxygen white dwarfs. The light curves and spectra we obtain 
from these simulations are in good agreement with observed properties of Type la supernovae. In particular, 
for white dwarf masses from 0.97-1.15 M Q we obtain 56 Ni masses between 0.3 and 0.8 M , sufficient to 
capture almost the complete range of Type la supernova brightnesses. Our optical light curve rise times, peak 
colours and decline timescales display trends which are generally consistent with observed characteristics al- 
though the range of B-band decline timescales displayed by our current set of models is somewhat too narrow. 
In agreement with observations, the maximum light spectra of the models show clear features associated with 
intermediate mass elements and reproduce the sense of the observed correlation between explosion luminosity 
and the ratio of the Si II lines at A6355 and A5972. We therefore suggest that sub-Chandrasekhar mass ex- 
plosions are a viable model for Type la supernovae for any binary evolution scenario leading to explosions in 
which the optical display is dominated by the material produced in a detonation of the primary white dwarf. 
Subject headings: radiative transfer — supernovae: general — white dwarfs 



1. INTRODUCTION 

In recent years, considerable work has been devoted to the 
study of the Chandrasekhar-mass (Mp h) model of Type l a 
supernovae (SNe la). As shown by lArnett et"a!1 (11971b . 
prompt detonations of Men carbon/oxygen (C+O) white 
dwarfs (WDs) in hydrostatic equilibrium mainly produce iron 
group elements (IGEs). Thus, they cannot account for the 
significant amounts of intermediate-mass elements (IMEs; 
e.g. silicon and sulphur) responsible for the features which 
dominate the maximum light spectra. To obtain these, pre- 
expansion of the WD is necessary such that burning partially 
takes place under low-density conditions where IMEs can be 
synthesized. One way of achieving this is provided by mod- 
els in which the flame ignites as a deflagration which releases 
sufficient energy to expand the star before a deflagration-to- 
detonation transition occurs (Khokhlov 1991). An alternative 
to this pre-expansion is the detonation of a sub-Chandrasekhar 
mass (sub-Mch) WD starting from a hydrostatic configuration. 
Here, a variety of density profiles can be realized, determined 
by the WD mass. Close to Ma,, the detonation produces pri- 
marily IGEs and few IMEs, while for less massive WDs more 
IMEs and less IGEs will be synthesized. 

Detonation of a sub-Ma WD cannot occur spontaneously 
but must be triggered by external compression. The most 
widely discussed mechanism for sub-Ma explosions has been 
the double detonation model. Here, a C+O WD accretes from 
a companion star and develops a helium-rich outer shell. This 
may occur for binaries with helium-rich donors or hydrogen- 
rich donors where the accreted hydrogen is burned to he- 
lium. If the helium-shell becomes sufficiently massive, it 
can become unstable and detonate. Subsequent compres- 
sion of the core by inward propagating shocks may pro- 
duce a secondary carbon detonation which exp lodes the WD 
(e.g. lWooslev & Weaver! 1986tlFink et aT]|2007l) . Detonations 



in helium-rich surface layers have also been discussed for 
the case of rapid dynamical mass transfer in binary systems 
containing a C+O WD with a helium-rich WD companion 
dGuillochon et al.f2010f) . In that case instabilities in the accre- 
tion seed dense knots which, by impacting on the underlying 
WD surface, might trigger a detonation in the accreted helium 
leading to a potential secondary core detonation. It has also 
been speculated that sub-Ma, explosions may arise during vi- 
olent accretion in mergers of C+O WD binaries. Here, the 
C+O accretion may lead to an edge-lit deto nation or carbon 
flashe s that trigger a core detonatio n (see e.g.lShigevama et alj 
1992 but for a different result see lLoren- Aguilar et al. 2009). 

Most previous work on testing sub-Ma models has fo- 
cused on cases in which the core detonation is triggered by 
detonation in an overlying massive shell (~ 0.2 M^) of he- 
lium (e.g. IWooslev & Weaverill994lLiyne& Arnettlll995l: 
Hoflic h & Khokhlovlll996t iHoflich et al.lll996t iNugent et alJ 
119971) . In those models burning in the helium shell syn- 
thesizes significant masses of 56 Ni in the outer ejecta, lead- 
ing to spectra and light curves in conflict with observations. 
As noted in those studies, however, these conclusions are 
strongly dependent on the influence of the shell material. In 
particular, they may not be applicable if such a layer is ab- 
sent (or much less massive) or if its post-burning composition 
lacks 56 Ni. Recently. [Bildsten et aTJ (120071) suggested that det- 
onation of the helium-shell may be possible for a shell with 
mass as low as ~ 0.055 M© around a 1 .025 M Q C+O core and 
that the burning pr oduces only 0.012 M© of 56 Ni along with 
some lighter IGEs (Guillo chon et al.l2010l find that even lower 
atomic-number burning products can dominate in their helium 
detona t ions). Even for th e low shell masses of Bildste net alj 
(2007), Fin k et all (|2010) find that a secondary core detona- 
tion is possible. To date, sub-Ma, explosions in the absence 
of a nickel-rich outer layer have not been studied in detail. 
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Shigeyama et al] ( I 19921) investigated the explosion dynamics 
of sub-Met, detonations and concluded that their characteris- 
tic properties were consistent with SNe la but they did not 
perform realistic radiative transfer simulations. 

A full treatment of any class of sub-Mo, explosion model 
requires realistic hydrodynamical and nucleosynthesis simu- 
lations of the accretion phase, triggering mechanism and sub- 
sequent explosion. Here, however, we present a simple nu- 
merical experiment that is relevant to any class of sub-Mo 
explosion model. We consider pure detonations of sub-Mo, 
C+O WDs with different masses, neglecting the question of 
how this detonation is initiated. This allows us to investigate 
the idealized case of sub-Met, detonation scenarios in which 
the observable display is dominated by material produced in 
the core explosion. Our goal is to determine the extent to 
which the least ambiguous component of the system, namely 
the detonation of a sub-Mch C+O WD, could lead to explo- 
sions which are consistent with observations of SNe la. 

2. HYDRODYNAMIC SIMULATIONS 

To explore the properties of detonations of sub-Met, WDs, 
we set up five hydrostatic models with WD masses (Mwd) 
ranging from 0.81 to 1.15 M© (see Tabled with centrally- 
ignited detonations. We performed the simulations in an ax- 
isymmetric setup of the full star but as the initial conditions 
were spherically symmetric and the symmetry is preserved 
during evolution, our models are one-dimensional. The sim- 
ulations were performed with our SNe la explosion code (the 
initial siz es of the compu tational grid cells are given in Ta- 
ble[2j see Fink et al. 2010 for technical details). 

The detonations were re presented with the level-set tech- 
nique (Reineck e et al.|[l999l) . This requires as inputs the deto- 
nation velocity and the energy release in the burning. The det- 
onation velocities take into ac count pathological d etonation 



speeds at high fuel densities (Gamezo et al.l 1 199*91) whereas 



low-density detonat ions are assum ed to be of Chapman- 
Jouguet type (see iFink et al. 2010 for details). The energy 
released by nuclear burning (E nuc ) and the asymptotic kinetic 
energy of the ejecta (E%) are given for each model in Table [2] 

3. NUCLEOSYNTHESIS CALCULATIONS 

The nucleosynthesis was computed using our standard 
tracer-particle technique (Travaglio et al. 2004). The masses 
obtained for 56 Ni, stable IGEs, IMEs and O are given in Ta- 
ble |2] 

We first computed nucleosynthesis for all five explosion 
simulations assuming an initial WD composition of uniformly 
mixed 12 C and 16 with equal mass fractions ("pure-C+O", 
hereafter). Although commonly adopted in Type la explo- 
sion modelling, this composition is not strictly correct since 
the mass-fraction of 16 Q is expected to b e larger than 12 C in 
the inner regions (e .g. ISalaris et alJI 19971). F or Met, delayed- 
detonation models, Dommgue z et al.l (1200 ll) showed that the 
C/O-ratio affects the 56 Ni-mass by ~ 14% and the veloc- 
ity structure of the ejecta by up to a few 1000 km s" 1 . Our 
calculations are expected to have a similar sensitivity to the 
adopted C/O-ratio. Moreover, C+O WDs formed from pro- 
genitors with non-zero metallicity will be polluted by some 
22 Ne. Since the neutron excess of 22 Ne significantly affects 
the nucleosynthesis, for one of the hydrodynamical mod- 
els (Mwd = 1.06 M©) we repeated the nucleosynthesis post- 
processing step adopting a high initial 22 Ne mass fraction of 
7.5% (this would corresponds to a rather high metallicity pro- 
genitor, Zq ~3Zq). This model ("C+O+Ne", hereafter) allows 
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FIG. 1 . — Composition versus mass-coordinate for the four models that pro- 
duced significant masses of 56 Ni (least- to most-massive; top to bottom). The 
black histograms indicate low-mass elements, blue shows intermediate-mass 
elements (Z = 1 1 — 20). Red represents 56 Ni while stable iron group elements 
(Z = 21—30) are shown in green. In the third panel, the solid lines show the 
results for the pure-C+O model while the dashed lines show the C+O+Ne 
model. In each panel the velocity scale (in Mm s ) is also indicated. 

us to bracket some of the systematic uncertainties associated 
with progenitor composition. The nucleosynthesis yields for 
this model are also given in Table [2] 

The range of Mwd we consider leads to 56 Ni-masses 
from ~ 0.01 M Q to 0.81 M Q , wide enough to encompass 
the range implied for all but the brightest SNe la (see e.g. 
IStritzinger et al.|[2006l) . The lowest mass models Mwd = 0.81 
and 0.88 M make very little 56 Ni (~0.01 and 0.07 M , re- 
spectively). Thus they would be faint and lie outside the range 
of normal SNe la. Throughout the following, we therefore ne- 
glect further discussion of the Mwd = 0.81 M Q model but re- 
tain the 0.88 M Q model as a point of reference for the faintest 
observed SNe la. 

Figure Q] shows the stratification of the nucleosynthesis 
products grouped into low-mass elements, IMEs, stable IGEs 
and 56 Ni. T he structure of the mod els is very similar to that 
obtained by Shigeya ma et al.l dl992l) : small masses of stable 
IGEs are produced and the mass-shell in which large mass- 
fractions of IMEs are produced is fairly extended. Significant 
IME mass-fractions are present up to almost the highest veloc- 
ities in all models, consistent with the lo wer limits on the oute r 
extent of Si-rich material discussed by Mazzali et al. (2007). 
Moreover, a clear trend exists whereby the inner boundary 
of the IME-rich layers lies at higher velocities in the models 
where the 56 Ni-mass is larger, the same trend as inferred from 
observations ( Mazzali et al. 2 007b . 

The most important consequence of 22 Ne in our C+O+Ne 
model is a substantial increase in the mass of stab le IGEs (see 
Table [2] and third panel of Figure [TJ see also Hofli ch et alj 
1998). These extend over a wide range of mass-coordinate 
and come at the expense of less 56 Ni in the inner regions and 
fewer IMEs in the outer zones. 

4. RADIATIVE TRANSFER SIMULATIONS 

For each of the four detonation models that produce > 
0.05 M Q of 56 Ni, we performed radiative transfer simulations 
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TABLE 1 

Nucleosynthesis products and optical light curve properties for detonations of white dwarfs. 



Parameter Models 



Kd < m 0) 


1.15 


1.06 


1.06 


0.97 


0.88 


0.81 


p c c (g crrT 3 ) 


7.9 x 10 7 


4.15 x 10 7 


4.15 x 10 7 


2.4 x 10 7 


1.45 X 10 7 


1.0 x 10 7 


WD comp. (C/0/Ne) rf 


50/50/0 


50/50/0 


42.5/50/7.5 


50/50/0 


50/50/0 


50/50/0 


(foe) 


1.64 


1.41 


1.41 


1.19 


0.97 


0.82 


E{ (foe) 


1.39 


1.22 


1.22 


1.04 


0.86 


0.73 


Cell size' (cm) 


1.0 x 10 6 


1.1 x 10 6 


1.1 x 10 6 


1.3 x 10 6 


1.5 x 10 6 


1.7 X 10 6 


Ml m (M ) 


0.81 


0.56 


0.43 


0.30 


0.07 


0.01 


Mf GE (Mq) 


2.1 x 10~ 2 


1.7 X 10~ 2 


1.8 x 10-' 


1.2 x 10- 2 


4.7 X 10~ 3 


8.7 x lfr 4 


Km ( M ©) 


0.27 


0.41 


0.36 


0.54 


0.63 


0.57 


M g Q (M ) 


0.04 


0.08 


0.09 


0.12 


0.17 


0.22 


Amis (mag.) 


1.34 


1.56 


1.42 


1.73 


1.77 




f Sax ( ma g ) 


19.2 


20.1 


18.0 


19.9 


14.1 




SLx ( ma g ) 


-19.9 


-19.2 


-18.7 


-18.5 


-16.6 




VL* (mag.) 


-19.6 


-19.4 


-19.3 


-18.8 


-17.3 




s Lx (mag.) 


-19.8 


-19.3 


-19.3 


-18.8 


-17.7 




'max (mag.) 


-19.6 


-19.0 


-19.2 


-18.7 


-17.8 




(B-V)Lx(mag.) 


0.15 


0.13 


0.48 


0.24 


0.63 




v Sffl (kms ') 


12,500 


11,500 


11,500 


9,000 


6,000 





" Since it has very low 56 Ni -mass, we did not perform radiative transfer simulations for the 0.81 Mq model. 

b Mass of white dwarf. 

c Central density of white dwarf. 

d Initial composition of WD (percentage by mass of I2 C/ lfi O/ 22 Ne) 

' Energy released by nuclear burning (£ n uc) and asymptotic kinetic energy of the ejecta (E^). 

' Initial size of computational grid cells in the WD. Since the simulations use an expanding grid, the physical 

resolution degrades with time. 

* Mass yields for 56 Ni, stable iron group elements (IGEs), intermediate-mass elements (IMEs) and oxygen (O). 
" Peak magnitudes are given at the true peaks in each band. Colours are quoted at time (f max ) of fi-band maximum. 
' Blueshift velocity of Si ii A6355 at t^. 



using our Monte Carlo code ARTIS dSiml2007 l :lKromer & Siml 
2009). For M WD = 1.06 M , we ran simulations for both our 
pure-C+O and C+O+Ne models. For all calculations we used 
our largest atomic data set ( ~ 8.2 x 10 6 lines) and our NLTE 
treatment of ionization (see iRromer & Siml 120091) . Table |2] 
gives the light curve decline-rate parameter (Amisfl the time 
of B-band maximum light (f^ ax ), the optical peak magnitudes 
and the B-V colour at f^ ax . 

5. COMPARISON WITH OBSERVATIONS 

In Figure [2] we show the ultraviolet-optical-infrared 
(UVOIR) bolometric and the band-limited (U, B, V, R, I, 
J, H and T^-band) light curves from our radiative trans- 
fer simulations. Maximum light spectra are shown for the 
same calculations in Figure [3] Observations of two SNe la 
(SN 2005cf and SN 2004eo) are shown for comparison in both 
figures. The 56 Ni-masses reported for these objects are signif- 
icantly different but withi n the range covered by our models 
(0.45 Mq for S N 2004eo llPastorello et al.l2007all and 0.7 M Q 
for SN 2005cf llPastorello et aUl2007bll ). 

Given the simplicity of the underlying explosion simula- 
tions, the light curve shapes and colours are in remarkably 
good agreement with observations. For Mwd = 0.97, 1.06 and 
1.15 M Q we obtain B-band rise tim es of 18-20 days, cl ose to 
observational estimates (~ 19 days; IConlev et alj|2006l) . The 
peak colours are also close to those observed but slightly red- 
der i n B-V compared to Mch models of si milar brightness 
(c.f. iHoflich et al.l [19961: iKaserTeTai] 12009). These results 
differ from previous studies fo r sub-Mch models where rel- 
atively rapid rise times (see e.g. Hoflich et al. 1996) and blue 

1 Amu is defined as the change in B-band magnitude between maximum 
light and 15 days thereafter. 



colours (c.f. Hofl ich et al.l 119961: iNugent et alj 119971) were 
found. Both systematic differences arise due to significant 
amounts of 56 Ni present in the outer layers of their models 
(see discussion by Hoflich & Khokhlov 1996!). Thus, our cal- 
culations illustrate that core detonations can be triggered 
without producing large masses of IGEs in the outer layers 
then good agreement with observations can be obtained. A 
modest additional mass of IMEs or unburned fuel ( 12 C, ie O 
or helium) in the outer ejecta due to the triggering mechanism 
would have minor consequences for the spectra. 

Our maximum light spectra are also in qualitatively good 
agreement with observations (Figure[3]). The models all show 
the characteristic Si II A6355 feature. Moreover, they repro- 
duce the sense of the observed trend whereby the strength 
of the weaker Si II feature at A5972 r elative to A6355 is 
systematically small e r in brighter events (|Nugent et alJI 1995b 
Bong ard et alJ l2006t lHachinger et al.l 12008). In agreement 
with observations, the maximum light spectra show clear fea- 
tures associated with other IMEs, in particular Ca and S. 
Our fainter models also predict O I absorption (A7773) as 
observed in some SNe la (including SN 2004eo; see Fig- 
ure [3]) but this feature becomes weak for our brighter mod- 
els, a trend also consistent with observations (Nuge nt et alJ 
119951) . As expected from Figure Q] there is a tendency for 
higher velocities of IME features in brighter events. Veloc- 
ities for Si II A6355 measured from our spectra are given in 
Tableland are gen erally compatible w ith those inferred from 
observations (e.g. iBenetti et alJ 1200 5). The Si velocity for 
our pure-C+O Mwd = 1 06 M Q model is slightly too high for 
both SN 2004eo and SN 2005cf. However, this discrepancy 
is small and on the scale of the differences between the mod- 
els (the observed line velocities are bracketed by those in our 
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FIG. 2.— Computed UVOIR bolometric, U, B, V, R, I, J, H and K-band light curves for M W d = 1.15, 1.06, 0.97, 0.88 M (red, green, blue and black, 
respectively). For Mwd = 1 -06 Mr , light curve s for both our pure-C+O and C+O+Ne mo dels are shown (solid and dashed green lines, respectively). Photometry 
for SN 2004eo (black circles; IPastorello etH] 20073) and SN 2005cf (black triangl es; IPastorello etal]|2007bl) are also shown. The observations have been 
corrected for reddening and distance using parameters from jPastorello et al. 20073lb|). 



Mwd = 0.97 and 1 .06 M Q maximum light spectra). 

Figure [4] shows the Z?-band width-luminosity relationship 
obtained from our models compar ed with the proper ties of 
a sample of well-observed SNe la (Hicken et al. 2009). The 
models reproduce the correct systematic trend: brighter mod- 
els have light curves which decline more slowly. Compared to 
the observations the more massive models (Mwd = 1-06 and 
1.15 Mq) have A/1115 larger than observed for their bright- 
ness. We note, however, that Am^ is a particularly chal- 
lenging quantity to model precisely and systematic uncertain- 
ties in the radiative transfer simulations can affect this quan- 
tity significantly. For example, applying diff erent radiative 
transfer codes to the well-kno wn W7 model ( Nomo to et al.l 
119841; iThielemann et al.ll 19861) . which is widely regarded as 
a good standard for normal SNe la, yields Am^-values that 
differ by several tent hs of a magnitude (see e.g. figure 7 of 
iRromer & Sim 2009). Moreover, the values obtained are also 
too large compared to those of normal SNe la (e.g. ARTIS 
yields A/M15 ~ 1.75). The predicted width-luminosity rela- 
tion is also affected by details of the initial WD as illustrated 
by our C+O+Ne test model. Compared to the equivalent pure- 
C+O model, this model declines more slowly in B-band and 
is fainter. This effect is sufficient to move this model to the 
opposite side of the observed width-luminosity relation (Fig- 
ure|4). Thus there may be potential for better agreement with 
observations from more detailed studies. 



Finally, we note that the near-infrared light curves ob- 
tained from the models have characteristic features which 
are in qualitative agreement with observations of SNe la. 
For the brighter models (M W d = 0.97, 1.06, 1.15 Mq), the 
/, J, H and K light curves show distinct secondary max- 
ima while the faintest model (Mwd = 0.88 M Q ) has single 
peaks. This is consistent with observations - normal SNe la 
show double-peaked near-infrared light curves while sub- 
lumin ous events have single maxima (e.g. IWood-Vasev et alj 
2008). Also, the difference between the models in the J, 
H and /^-bands at B-band maximum light is much smaller 
than in the optical bands, consistent with the observation 
that SN e la are better standard candles at near-infrared wave- 
lengths dKrisciunas et al.ll2004 . 

6. DISCUSSION 

The sub-Mch model for SNe la has much to commend it. 
First, it has already been suggeste d by empirical mode lling 
of bolometric SNe la light curves (Stritzinger et al. 2006)) that 
differing ejecta masses may be required for different SNe la, 
a property which the sub-Mch model may explain. Secondly, 
population synthesis studies predict large numbers of binary 
systems with accreting C+O WDs: iRuiter et ail (120091) esti- 
mate a Galactic rate of ~ 10~ 3 yr" 1 for possible explosions of 
sub-Mch C+O WDs accreting from helium-rich companions. 
This is much higher than their estimate of the Galactic rate 
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FIG. 3 . — Computed spectra at B-band maximum light for Mwd = 1.15, 1 .06, 0.97, 0.88 Mq (red, green, blue and black, respectively). For Mwd = 1 06 Mq , 
spectra for both our pure-C+O and C+O+Ne models are shown (thick an d thin green lines, respectively). Observ ed spectra of two SNe la around maximum light 
are shown for comparison: SN 2004eo (lower grey line; Pastorello et al. 2007a) and SN 2005cf (upper grey line; Garavini et al. 2007). Arbitrary vertical offsets 
have been applied for clarity. The observed spectra are de-redshifted and de-reddened using parameters from Pastorello et al. l2007ajJ3). 

ample, if it can be shown that explosions in binary systems 
with larger Mwd are more often found among young stellar 
populations relative to their less massive Mwd counterparts, 
the observ ed correlation of SN la brightness with host galaxy 
type (e.g., Howell 2001) might be explained. 

Here we have shown that detonations of sub-Met, WDs lead 
to explosions which give a reasonable match to several prop- 
erties of SNe la. Specifically, WDs with masses between ~ 1 
and ~ 1 .2 M Q can reproduce a wide range of brightness with 
light curves that have rise times and peak colours in roughly 
the correct range. In addition, the models reproduce the char- 
acteristic spectral features present around maximum light and 
the observed trend for a higher velocity at the inner bound- 
ary o f the IME-rich layer in brighter SNe la dMazzali et al.l 

2007) . Although our pure-C+O models yield light curves 
that fade too fast after maximum, the models predict a width- 
luminosity relation which behaves in the observed sense and 
we would argue that the combination of uncertainties in ra- 
diative transfer simulations and details of the nucleosynthesis 
(which is sensitive to the progenitor composition) can system- 
atically affect the decline timescale. Thus there is potential for 
even better agreement with improved modelling. 

There are several additional observational constraints that 
our current models do not address but which should be con- 
sidered in future studies. For example, off-centre detonation 
might lead to observable effe cts associated w ith departures 
from spherical symmetry (e.g. Fin k et alJ2010l) . Chemical in- 
homogeneity of the pre-explosion WD could affect the explo- 
sive nucleosynthesis: in parti cular, sign i ficant gravitational 
settlin g of 22 Ne dBildsten & Halll I200U iGarcfa-Berro et alJ 

2008) might yield a layered ejecta structure with a central 
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FIG. 4. — Comparison of the B-band light curve peak magnitude (B max ) and 
decline-rate parameter (Amu) rel ation obtained from the models (solid cir- 
cles) with observations of SNe la iHicken et al. 2009; grey crosses). SNe la 
with distance modulus /i < 33 have been excluded. The open circle shows 
the C+O+Ne model and is connected by a dashed line to the point for the 
pure-C+O model with the same Mwd ■ 



for single-degenerate Met, explosions (0.6- 1.4 x 10 yr -1 ) 
and comparable to their estimate of the WD-WD merger rate 
(1 - 2 x 10~ 3 yr" 1 ) in systems that exceed Men- For compari- 
son, the observed Gal actic rate of SNe la is (4 ± 2) x 10~ 3 yr -1 
(ICappellaro et al.ll 19991) . 

Moreover, sub-Mch models provide a simple physical pa- 
rameter which could account for the range of observed bright- 
nesses: the mass of the exploding C+O WD. This parameter 
allows for a possible link between the typical brightness of a 
SN la and the stellar population in which it resides. For ex- 
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concentratio n of neutron-rich isotopes as favoured b y obser- 
vations (e.g. iHoflich et al.ll2004tlGerardv et alJl2007l) . 

In conclusion, detonations of naked sub-Met, C+O WDs 
yield light curves and spectra which are in qualitatively good 
agreement with the observed properties of SNe la. The critical 
question remains whether or not realistic progenitor scenarios 
in which the optical display is dominated by such an explosion 
can be established: it must involve detonation of a WD with 
a density profile similar to those of our toy models without 
producing large masses of high-velocity IGEs. Any sub-Mch 
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